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Abstract 

In the framework of Left-Right symmetric model, we investigate an interesting scenario, in which 
the so-called VEV seesaw problem can be naturally solved with Z2 symmetry. In such a scenario, 
we find a pair of stable weakly interacting massive particles (WIMPs), which may be the cold dark 
matter candidates. However, the WIMP-nucleon cross section is 3-5 orders of magnitude above 
the present upper bounds from the direct dark matter detection experiments for m 10^ - 10'* 
GeV. As a result, the relic number density of two stable particles has to be strongly suppressed 
to a very small level. Nevertheless, our analysis shows that this scenario can't provide very large 
annihilation cross sections so as to give the desired relic abundance except for the resonance case. 
Only for the case if the rotation curves of disk galaxies are explained by the Modified Newtonian 
Dynamics (MOND), the stable WIMPs could be as the candidates of cold dark matter. 
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I. INTRODUCTION 

The Left-Right (LR) symmetric model [ll, based on the gauge group SU{2)l x SU{2)r x 
U{1)b-Li is an attractive extension of the standard model (SM). The symmetry requires 
the introduction of right-handed partners for the observed gauge bosons and neutrinos, and 
a Higgs sector containing one bi-doublet (2,2,0), one left-handed triplet A/, (3,1,2) and 
one right-handed triplet A/j (1,3,2). In such a minimal LR symmetric model, parity is an 
exact symmetry of the theory at high energy scale, and is broken spontaneously at low 
energy scale due to the asymmetric vacuum. Also CP asymmetry can be realized as a con- 
sequence of spontaneous symmetry breaking, namely the spontaneous CP violation (SCPV) 
{2!. However, such a scenario suffers from nontrivial constraints from the vacuum minimiza- 
tion conditions. It is explicitly demonstrated that the SCPV is not so easily realized if all 
the parameters in the Higgs potential are real and endowed with natural values 3|, |J, |5|. The 
difficulty results from the facts that one of the neutral Higgs bosons carries dangerous tree 
level flavor changing neutral currents (FCNC) effect, and that quark flavor mixing angles 
and CP violating phase are all calculable quantities due to the LR symmetry. Therefore, 
many generalized CP violation scenarios beyond the SCPV case have been analyzed ex- 
tensively , [lOl, lUl- In these literatures, the masses of right-handed gauge boson 



W2 and the FCNC Higgs boson are strongly constrained from low energy phenomenology. 
Although the CKM matrix are more general not to be fully fixed than the SCPV case, it is 
proved that there is only one physical complex phase in the Yukawa couplings [l2[. Hence 
the FCNC Higgs boson's couplings can't be absolutely free. The FCNC Higgs boson's mass 
still accepts strict bound. In terms of these observations, a generalized two Higgs bi-doublets 
model is proposed [l^. In this model, quark mass matrices become far more fiexible and 
the FCNC Higgs boson's Yukawa couplings are now free parameters. Thereby low energy 
bound on the right-handed scale is largely alleviated. As other generalized models, the two 
Higgs bi-doublets version of LR model also has the advantage to realize the SCPV without 
the fine-tuning problem. 

The LR symmetric model is also motivated to explain the very tiny neutrino masses. 
When the vacuum expectation value (VEV) Vr of the neutral component of A/? is very 
huge, typically of order 10^^ GeV, the well-known seesaw mechanism provides a very natural 
explanation of the smallness of neutrino masses [l^. However, the right-handed gauge 
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bosons Z2 and Wi are too heavy to be detected at the Large Hadron Colhder (LHC) and 
the future coUiders. To allow for the possibility of an observable right-handed scale, many 
authors focus on the f/j ~ 10 TeV case. Although the seesaw mechanism can work well, we 
have to face the so-called VEV-seesaw puzzle. Namely, is of order 10^^*^ rather than 
the anticipant 0{\\ where p and (5 are located in the Higgs potential. One may introduce 
a discrete Z2 symmetry A/, — and A/? Ar to resolve this VEV-seesaw problem 
^. It is worthwhile to stress that neutrinos are the Dirac particles in this scenario. If we 
preserve the Majorana Yukawa couplings, the corresponding model must lie beyond the LR 
symmetric model. 

The Z2 symmetry leads to the absence of both /3-type terms and the Majorana Yukawa 
couplings, hence f ^ = due to the minimization conditions. Furthermore, we find that the 
neutral Higgs bosons b\ and 5^* are a pair of stable weakly interacting massive particles 
(WIMPs). This is an important feature of our scenario which hasn't been indicated before. 
It is a natural idea that 6\ and 5^* may be the cold dark matter candidates |l5|. We 
firstly calculate the WIMP-nucleon elastic scattering cross section which has been strongly 
constrained by the direct dark matter detection experiments, such as the CDMs[l6| and 



XEN0N[17l]. However, our result is 3-5 orders of magnitude above the present bounds for 
m ~ 10^ - 10^ GeV 
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171 ]. To avoid this puzzle, 6\ and can't dominate all the 



dark matter. We find that our scenario is consistent with the direct dark matter detection 
experiments only when ra^o < 4.8 x 10"^^, where n^o is the total relic number density of 
b\ and This bound requires the dark matter annihilation cross sections must be very 
large. In this work, we examine whether our scenario can provide very large annihilation 
cross sections so as to derive the desired relic abundance. 

In this paper we try to give a comprehensive analysis on these LR models with general 
parameter setting. Firstly, we perform a detailed investigation on the simplest LR model 
with one Higgs bi-doublet, in which there are no any CP violation phases. Then we generalize 
the simplest LR model to some other more complicated situations. It turns out that there's 
no significant differences among these one Higgs bi-doublet versions of LR model because 
the gauge and Higgs sectors are basically the same. Whereas in the two Higgs bi-doublet 
case, there would be more Higgs bosons and the Yukawa couplings might be quite different. 
Hence more delicate analysis is needed. The remaining part of this paper is organized as 
follows. In Section II, we briefly describe the main features of the LR symmetric model and 



3 



discuss the VEV-seesaw problem. In Section III and IV, the direct dark matter detection 
experiments put very strong constraints on the relic number density and the annihilation 
cross sections. In Section V, we analyze whether the simplest LR model can be consistent 
with the above constraints or not. Then we generalize the simplest LR model to the two 
Higgs bi-doublets case in Sec VI. The summary and comments are given in Section VII. 



II. THE LR SYMMETRIC MODEL WITH SYMMETRY 



The minimal LR symmetric model consists of one Higgs bi-doublet (2,2,0), one left- 
handed Higgs triplet (3,1,2) and one right-handed Higgs triplet A^ (1,3,2), which can 
be written as 

After the spontaneous symmetry breaking, the Higgs multiplets can have the following 
vacuum expectation values 

\ / 



(2) 



where ki, K2, vl and vr are in general complex. Without loss of generality, one can choose 
Ki and vr to be real, while assign complex phases 62 and Ol for k2 and vl, respectively. 
Following the requirements of the LR symmetry, we can write down the most general form 



of the Higgs potential 
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where = T20*r2 and all parameters /ij, Aj, p,, and /3j are real. Only 02 can be complex. 
The phases of K2 and vl may lead to the SCPV It has been shown that the combing 
constraints from K and B system actually exclude the minimal LR symmetric Model with 
the SCPV in the decoupling limit j^. For our present purpose, we investigate here the 
simplest LR model, in which 02, ^2, and the Yukawa couplings are real. It is worthwhile 
;o stress that our remaining analysis can be generalized to the other CP violation scenarios 
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In the minimal LR symmetric model, the Lagrangian relevant for the neutrino masses 
reads 



(4) 



where ipL,R = {^l,rJl,r)'^ ■ After the spontaneous symmetry breaking, one may obtain the 
effective (light and left-handed) neutrino mass matrix m,^ via the type II seesaw mechanism: 
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where k = y^l^iP + 1^2^ ~ 246 GeV represents the electroweak symmetry breaking (EWSB) 
scale and Yd = (Y^jKi + YjjK2) / {\/2k). The charged lepton mass matrix is given by nii = 
{YyK2 + YyKi) / \/2. The electroweak precision test requires Vl <^ k. Barring extreme fine- 
tuning, the neutrino masses nii, ~ 0.1 eV 13] forces vl to be of order a few eV or less, 
thereby requiring vr ~ 10^^ GeV for Yd ~ Ym ^ mi/ k. In this case, the right-handed gauge 
bosons Z2 and W2 are too heavy to be detected at the LHC and the future colliders. To allow 
for the possibility of an observable right-handed scale, many authors focus on the vr ~ 10 
TeV case. Although the seesaw mechanism can work well, we need to resolve the so-called 
VEV-seesaw puzzle j4|, which is indicated by a simple vacuum minimization equation: 



(2pi - P3)vlVr = PiniK2 + p2i^l + P3l^l ■ 

Without loss of generality, one can write Eq.(l6]) in a compact form: 

P VlVr 



(6) 
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In view of the naturalness, one expects 7 ~ However, we find that 7 ~ 

long as ~ 10 TeV. This is the infamous VEV-seesaw problem in the literatures 
neutrino mass matrix rrii, in Eq. can also be written as 

.2 



O-^'J as 
The 
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It is shown that the VEV-seesaw relationship implies the unnaturalness for the auxiliary 
parameter 7 if one wants to search for new physics at TeV scale. To avoid the VEV-seesaw 
puzzle, a smart way is to introduce some new symmetries to eliminate all /3-type terms of 
the Higgs potential. However this is not a easy task in the current model. One may guess 
there exists some additional global symmetries like U{1) acting on the Higgs fields which 
can eliminate all /3-type terms ^. However, such alternative always affects the fermion 
sector and fails to give correct fermion masses and mixing. If there is an approximate U{1) 
horizontal symmetry to suppress 0^ without eliminating them completely, then one may 



solve the VEV-seesaw problem 10, ll9[. Unfortunately, this model yields a small mixing 
angle within the first two lepton generations. In Ref . J] , the authors suggest a Z2 symmetry 

Al ^ -Ai, Ar^Ar, (9) 

which can eliminate all /3-type terms of the Higgs potential. However, this discrete symmetry 
also eliminates the Majorana Yukawa couplings, which implies that neutrinos are Dirac 
particles. At this moment, Eq.Q becomes 

(2pi - P3)VL = . (10) 

One may immediately dismiss the possibility 2pi — = 0, which implies two massless 
left-handed Higgs triplet bosons. Thus the only left choice is vl = 0. The Z2 symmetry 
leads to Vl = and the absence of both /3-type terms and Majorana Yukawa couplings. 
Furthermore, we find that the lightest particles among the members of left-handed Higgs 
triplet Al, namely 6l and 6l*, are two degenerate and stable particles. A natural idea 
is that 6l and 6l* may be the cold dark matter candidates. In the following sections 
we shall discuss the possibility of 5^ and 5°* being the cold dark matter candidates by 
evaluating all relevant annihilation processes. The main features of the LR symmetric model 
with Z2 symmetry have been shown in Ref. [2^. Here, we show the mass spectrum for 
the Higgs bosons and gauged bosons at leading order in Table. I, with approximations 
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TABLE I: The mass spectrum for the Higgs bosons and the gauged bosons in the LR symmetric 
model with symmetry. Here, we have neglected the terms in order of ^2/^1 and /v\. 

n'^/vj^ ~ and 1^2/1^1 — mentioned in Appendix A. Gauge bosons Zi and Z2 are defined 
by Zi = cwWsL — swtw^sR — \/C2wtwB and Z2 = ^Jc2w secvy W^^r — t^B, where the 
subscript W denotes the Weinberg angle 9w- In addition, all the trilinear and quartic scalar 
interactions and scalar-gauge interactions are listed in Appendix A for convenience. 



III. THE DIRECT DARK MATTER DETECTION 



The current direct dark matter detection experiments, such as the CDMS[16| and 



XENON 



171 ]. have provided very strong constraints on the WIMP-nucleus e 



tion. The rate for direct detection of dark matter candidates is given by 15 1 



astic cross sec- 



Plocal 



m 



(11) 



where iVj is the number of nuclei with species i in the detector, piocai is the local energy 
density of dark matter, m is the mass of cold dark matter. aij\f is the WIMP-nucleus elastic 
cross section, and the angular brackets denote an average over the relative WIMP velocity 
with respect to the detector. Using the standard assumptions of piocai and distribution of the 



relative WIMP velocity 2l|, one can derive the constrains on WIMP-nucleon cross-section 
< 4.6 X 10-^^cm2 for m = 60 GeV from the CDMS Q; afJ' < 8.8 x 10"^^ cm^ for 



m 



100 GeV from the XENON 



17| . Since the WIMP flux decreases oc 1/m, o"^^'^ oc m is a 



very good assumption for m > 100 GeV. 



In our scenario, the dark matter candidates 5l and 5l* interact with nucleus J\f through 
their couphngs with quarks by exchanging the neutral gauge bosons Zi, Z2 and Higgs bosons. 
We find that the main contribution comes from the Zi exchanging process, which produces 
a spin-independent elastic cross section on a nucleus Af [22 1 

2GlM\Af). 
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where Z and A — Z are the numbers of protons and neutrons in the nucleus, respectively. 

is Fermi coupling constant and M{M) = mMj^/{m + Mj^) is the reduced WIMP mass. 
Traditionally, the results of WIMP-nucleus elastic experiments are presented in the form of 
a normalized the WIMP-nucleon cross section cr„ in spin-independent case, which is straight 
forward 

1 M'^in) 



(Jr. 
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M^{U) 

where M{n) = mMn/{m + M„) and M„ denotes the nucleon mass. When m 3> M„, one 
may arrive at cr„ = 8.2 x lQ~^^cm^ for the CDMS experiment^jwhich is 3 — 5 orders of 

10^ - 10^ GeV |l7|. Therefore, such dark 



magnitude above the present bounds for m ~ 10^ — 10"^ GeV 
matter candidates are excluded by the current direct detection experiments. 

If 5^ and have a nonzero splitting, one can avoid the above bounds since the Zi 
exchanging process is forbidden kinematically 23|]. However, such degeneracy can not be 
satisfied in our model. If the energy density of 5° and 5^^^ in the solar system is far less than 
Piocaii we can avoid the above experimental limits as shown in Eg. (ITT]) . This means that 
5\ and 5^* are only a very small part of the total dark matter. We find that our model is 
consistent with the direct detection experiments only when 



n^o < 4.8 X 10 
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where n^o is the total relic number density of 5\ and 5^*. Here we have taken the approx- 

1 TeV) as 



imation a'^^ oc m (when m > 100 GeV) and used o",^^^ = 3.4 x 10 cm 



m 



3- 



It is worthwhile to stress that the bound in Eg. 0141) is not valid 



the input parameter 
for m < 100 GeV. 

The present experimental bounds are based on the standard assumptions for the galatic 
halo 2l|. It needs to be mentioned that the rotation curves of disk galaxies may also be 



explained by the Modified Newtonian Dynamics (MOND) 2J]. On one hand, we use the 



MOND to account for the rotation curve of the Milk Way; On the other hand, we still believe 
that the cold dark matter exists in the universe. In this case, the local energy density of 
cold dark matter may be far less than the standard assumption. Therefore, we may give up 
the above constraints from the direct dark matter detection experiments. Subsequently, the 
stable particles 5° and 61* may be the cold dark matter. 



IV. CONSTRAINTS ON THE ANNIHILATION CROSS SECTION 



The thermal average of annihilation cross section times the "relative velocity" (av) is 
a key quantity in the determination of the cosmic relic abundances of 6^ and 6^* . The 
constraint in Eq. flT^ implies (av) must be very large in our scenario. In this section, we 
analyze whether the present model can satisfy Eq. lfT^ . 

In our scenario, S\ (i = 1, 6 for 6l,6l*,6t and 5f ^) set of similar particles whose 
masses may be nearly degenerate. The total relic density of the lightest particles 6^ and S^* 
is determined not only by their annihilation cross sections, but also by the annihilation of 
the heavier particles, which will later decay into 5° or Therefore, we need to consider 



the coannihilation processes [25|. Since 5^ and which survive annihilation eventually 



decay into 5° or 6^* , the relevant quantity is the total number density of 6\^, n = J2 



25|: 



i=l 



Hi. 



The evolution of n is given by the following Boltzmann equation 

^ = -2,Hn-{a,ffv){n'-n%), (15) 

where H is the Hubble parameter, Ueq is the total equilibrium number density, v is the 
relative velocity of two annihilation particles. The effective annihilation cross section cTg// is 

^e// = E^..^(l + ^^f'\l + A,)3/2e-^(^'+^^) , (16) 

where Aj = (mj — m)/m, x = m/T is the scaled inverse temperature. Qi = 1 is the internal 
degrees of freedom of 5]^ and geff = Z)f=i 5'i(l+Aj)'^/^e~^^\ For the total equilibrium number 
density, we may use the nonrelativistic approximation Ueq ~ (7e//(mT/27r)^/^exp(— m/T). 

For particles which potentially play the role of cold dark matter, the relevant freeze-out 
temperature is Xf = m/T ~ 25. In our scenario, one can derive Xf ^ 35 which can be seen 
in Eq. (fT9!l . When Aj > 0.1 for Sf; and 5^^, we can arrive at cTg/-/ = 0"i2/2 in our model. 
In addition, we find that it is also a rational approximation a^jf ~ o'i2/2 even if all masses 



of 5\ are nearly degenerate. For simplicity, we take ((Je//f) = {auv) /2 in the remaining 
analysis of our paper. 

For nonrelativistic gases, the thermally averaged annihilation cross section {cri2v) may be 
expanded in powers of (cuv) = (TqX^^ , A; = for the s-wave annihilation and k = 1 for 



the p-wave annihilation 26|]. The general formula for {(Ti2v) is given by [27|] 



(cri2t') =(ToX = ^ 



3 



, ' (17) 



where u = EiE2(Ji2V, prime denotes derivative with respect to s/4m^, and s is the center- 
of-mass squared energy, u and its derivative are all to be evaluated at s/Am^ = 1. The final 



number density n^o is given by 26 1 



3.79(A; + 1)4^^ 



ngo = 2970 ' ^ cm"^ (18) 



with 



Xf = \n[Om8{k + l){geff/9l/^)Mpimao/2] 

- (k + l/2)ln{\n[0m8ik + l)ig,ff/gl/^)Mpimao/2]} , (19) 

where Mpi = 1.22 x 10^^ GeV and is the total number of effectively relativistic degrees 
of freedom at the time of freeze-out. Here we take g^, ~ 100 for illustration. With the help 
of Eqs. ([HD, ([HD and ([ISD, we can derive 

mcTQ ^ 0.13 GeV^^ (s— wave) ; 

m(To > 9.8 GeV"^ (p-wave) . (20) 

V. ONE HIGGS BI-DOUBLET MODEL 

In this section, we shall investigate whether the above bounds can be satisfied in one 
Higgs bi-doublet model or not. Since there are many unknown parameters, some rational 
assumptions have to be made for our model so that one can calculate all relevant annihilation 
processes. In our scenario, the thermally averaged annihilation cross section {(Tuv) is usually 
inverse proportional to as shown in Eq. ( fT71) . Therefore, one can obtain mao oc I/vr. 
Namely, the smaller vp is, the easier Eq. (!20|) can be satisfied. Considering the constraints 



on the masses of W2 and the FCNC Higgs boson from low energy phenomenology 
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ll|, we 



choose vr = 10 TeV and as = 2 as an instructive example to illustrate the main features 
of our scenario. One can immediately get = 7.5 TeV and = 4.5 TeV. Now let's 
introduce an auxiliary parameter e = {ps — 2pi)/(2pi) to reexpress the mass of 5° and 5°* 



1^= ]j^{p3- 2pi) Vr = y/eplVR . (2l) 

From the Zi invisible width one may obtain m > mz^/l, which requires epi > 2.0 x 10"'^. 
On the other hand, we may require < 4 in view of the perturbativity, and then derive 
Pi + ^Pi ^ 2. In addition, we wish all pj have the same order which means s < 4. Due 
to the suppression of phase space, one may ignore some annihilation processes in terms of 
the values of e and pi. When m < mwi, 5£ and 5°* mainly annihilate into the fermion 
pairs (except for top quark). The corresponding mo"o is far less than the lower bound of 
Eq. fl20l) . For the convenience of the remaining analysis, we require m > 500 GeV (Namely, 
^Pi > 2.0 X 10"'^) which does not affect our conclusions. Finally, we assume that all ctj of 
the Higgs potential have the same order. 

It is worthwhile to stress that Eq. ( |T71) is not valid when the annihilation takes place near 



a pole in the cross section 25[. This happens, for example, in Z-exchange annihilation when 
the mass of relic particle is near m^/2. For the cases 2m/mz < 0.8 and 2m/mz > 1.2, we 
use the above analytic way to calculate mao. On the contrary, we should numerically solve 
the Boltzmann equation in Eq. ffTSj) . in which the resonant cross sections of the Breit-Wigner 
form must be considered. Then one can derive the relic number density n^o which has to be 
less than the upper bound in Eq.f ll4l) . 

In general, all relevant annihilation processes may be divided into four categories in terms 
of the different final states: S^Sl* ff, 5l5l* VV , HH and VH, 

where V and H denote the gauge boson and the Higgs boson, respectively. Next, we shall 
analyze in detail the four classes of annihilation processes and the resonance case. 



A. 5l6l*^ff 

Let's start with the first case: 5° and 6^* annihilate into fermion pairs. There are two 
kinds of Feynman diagrams at the tree level contributing to this case: S channel gauge bosons 
exchanging and Higgs bosons exchanging diagrams. Because of the absence of Majorana-type 
Yukawa couplings, there are no the T channel diagrams' contribution. The first amplitude 
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is proportional to e^, while the second is proportional to amj/^/s. It is plausible that both 
diagrams have the same contribution for m ~ 500 GeV. However, the squared amplitude 
of the first diagram always includes a suppression factor of 1 — Am? / s, which leads to the 
p-wave annihilation. 

For the gauge bosons exchanging diagram, we can obtain 

It is obvious that this is a p-wave annihilation process. With the help of Eq. (fT7|) . we have 
(mcTo)// < 2.2 X 10~^ GeV~^ for m > 500 GeV. It is 6 orders less than the lower bound 
> 9.8 GeV~\ Ahhough one may increase mao through lowering m, {mao)ff is still far 
less than the lower bound in Eq. fl^U]) even if m = 100 GeV. Therefore, this process can not 
suppress the relic number density of 5° and 5°*. 

For the Higgs bosons exchanging diagram, the exchanged particles should be /i° and Hi. 
As shown in Table I, the mass of is far more than the light SM Higgs mass rriho. Due 
to the suppression of propagator, we neglect the contribution from H^. For the case, the 
amplitude of Higgs bosons exchanging process is proportional to mj. Furthermore, we only 
consider the top quark pair final states. The relevant cross section is 

_ 3 {aimtYs . 

which leads to a s-wave annihilation process. One may immediately derive (mo"o)top ^ 
1.5 X 10"^ GeV~^ for m > 500 GeV and ai = 2, which is far less than the lower bound 
mao > 0.13 GeV~\ 



B. 6l6l* VV 

In Fig. [H we show all possible Feynman diagrams for the process 6^6^* — ^ VV. There 
are three kinds of Feynman diagrams Fig. la, lb and Ic for the final states ZiZi. Obviously, 
the amplitude of Fig. lb is suppressed by a factor of i^l\fs compared with the first one. 
Thus we only consider the contribution from Fig. la and Ic. The total annihilation cross 
section is found to be 

2e^ csc^ 2ew 

TT 

12 



1 + 



Am^ 8 



m^is 



2m^ 



(24) 



{b) 



{d) 



FIG. 1: All possible Feynman diagrams for the annihilation processes 6^6^* —>■ VV, where 
may be S^^^* or 6f^, and denotes /i*^, and 



where the function y{xi) is defined by y{xi) = arctanh(xi)/xi and xi = yl — Am? / s. Then, 
we can derive {maQ)ziZi ^ 2.1 x 10"^ GeV~^ for m > 500 GeV. It is obvious that this result 
is not so large as to satisfy the requirement of Eq. fl20|) . 

According to the ZiZi experience, we also calculate the other processes. The correspond- 
ing cross sections are given by 
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(25) 



(26) 



(27) 



(28) 



where Xp = (cot^ OwPsvji) / (s — m^o) and X2 = \J {s — Am'^)(s — Am'zJ/{s — 2m\^). These 
cross sections have the same order as the ZiZi case. However, the thermally averaged 
annihilation cross sections of these processes (except for Wiiyi) are far less than the ZiZi 
case with m = 500 GeV. Therefore, we don't analyze these processes in detail. 
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C. 6lSl* HH/VH 



X 



HO 



/ \ 



/ 



(a) 



Y 

Axj [T « U) 



FIG. 2: All possible Feynman diagrams for the annihilation processes S^S^* — > HH. 




FIG. 3: All possible Feynman diagrams are shown for the annihilation processes S^S^* — > VH. 
The first diagram only appears in the process 6^5 — > VA^ 

Let's now focus on the processes S^Sl* — > HH/VH. The relevant Feynman diagrams for 
HH and HV are shown in Fig. [2] and Fig. [3l respectively. Since the dimensional scalar tri- 
linear couplings enter extensively into the above two annihilation processes, the electroweak 
scale coupling ain in 6l6l*h^ and the right-handed scale coupling PsVr in 6i^6i^*H2 would 
make big difference in the 5° 5£* — > HV processes according to our current parameter setting. 
Considering the complexity of this model, we only calculate the annihilation cross sections 
up to leading order (LO) by omitting the next to leading order (NLO) contributions in 
terms of the following three suppressing factors: (1) small VEV ratio n/vfi and ^^/^/s due 
to the big hierarchy in symmetry breaking scale of the LR model. Since we have made the 
approximation n"^ /v\ ~ thus here it is of course a reasonable power counting rule to pick 
out the LO processes against the NLO ones; (2) gauge coupling suppression e^; (3) p-wave 
factor 1 — Am? / s due to large suppression in the integration of initial energy of dark matter 
pair. 

In this subsection, we apply the above three suppressing factors to make an explicit 
demonstration of the LO processes, then give the convincing dark matter annihilation cross 
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Process 



Amplitude 



Order 



Process 



Amplitude 



Order 



Him 



iai 1 



etw [c2wPzi - twPz2/'A 



ip3 2 



-2ip^a2KVR 



t—m^ 



1 J 1 J 1 

s—2piv'-L t—m'-^ u~m-^ 



iP3 



1 J 1 

s—2piv'i, t—m'^ 



^ - + ^^s-2pi4 



+ 



1 



2(pi+2p2H 
s~2piv'j^ 



P4 ' 



P3 t-rri'^ 



' K \2 



{ — ? 



H'lZ^ 

A\Z^ 
A\Z2 

2°^i 

2^2 



HnZo 



Aieai esc 29\ykPtu 
8iea2 csc2vk t^Pru 



Aiea2 csc2VK i^PE, 
-Aiea2 csc2Ty ^C2wi^PE 
4iep3 CSC2VF vrPtu 

—2iea2 csc2vy ^P^ 
—2iea2 csc2vy i^P^s 



Xf s 
Xf s 



(-)' 



e 



(-)' 



TABLE II: The amplitude for HH/HV final states, where cw = cos^iyj tw = tan^iy, etc. We 
also estimate the order of corresponding annihilation cross sections. 



sections. The LO amplitude for each possible annihilation process is listed in Table |TT1 The 
notations are as follows: pi^2 denotes the momentum of the dark matter pair, while ^3^4 is 
the momentum of the final states, e is the polar vector of gauge boson: 



(Pi - V2) ■ {Va - Ps) 



n,2 



m 



-jh,H 
43 



Pa ■ ^{Ps) 



^1,2 



m 



p P2 ■ e(P3) Pi ■ e{p3) 

ptu - — — — \ — — i-^yj 



hO,H<l 



t — 



u — 



Here we only consider the cross sections with amplitude order 1. In terms of Table [TTl 
nine LO annihilation cross sections are listed in Table IIIIl where A = 1 — ^(pi+^p^) 
v(x-i) + — 2 ^ 2 — ; — 3 — with X'i 

3m? „ 



B 



pi s—2m^,. 



„ r- and 

s-2p\vj^ 



s — Am?){s — 4m^±±)/ (s — 



4m^ + miro , c 
"2 



s — 2m?,o, d 
"2 



m^ro + rn^jo and 



X4 = ^{s — Am?){s — 4m^o)/ (s — 2m^o). We find that these processes fail to provide enough 
large cross sections. For 52, 5^,* — > A\Zi and HfWi, one can easily obtain mao ^ a\/ (167rm), 
which is far less than the required lower bound 0.13 GeV~^. Since the other processes have 
the similar forms, we take the process ^'j^d^j^ h^H^ as an example to illustrate the main 
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Process 



AEiE2(JV 
#(1 



Am'' 



P3ff 



s—2p\vi 



Am'' 



167r 



(1 



s—2p\vi 



167r 



a2 + ^ + f(2a + ^Mx4) 



Process 



AEiE2(Tv 



+ ^ 



oi ' s—2p\v'i 



4m' 



±± 



1/2(^2^5) 



.)5/2 



)5/2 



TABLE III: The annihilation cross sections for the leading order processes, 
features of this kind of processes. One can immediately derive 



(mao) 



1 



1 



Sttvr^ epi A{epi)^ \2e-lJ ' ^^^^ 
where we have used as = 2. It is obvious that the maximum value can be obtained when 
epi = 0.5. Varying e, we may derive {'mao)f^oHO < 3.5 x 10"^ GeV""*^ for a2 = 2. Therefore, 
we don't discuss this class of processes in detail. 



D. The resonance case 



As pointed out in the previous discussion, the method of calculating the effective t 
mally averaged annihilation cross section {ceffv) is not valid for the resonance case 



ner- 



25| 



3ere we numerically solve the Boltzmann equation Eq. ffTSj) . which can be reexpressed as 



28| 



dY 

dx 



X 



y2 



(31) 



where H^^i is the Hubble parameter evaluated at T = m and s(x) is the entropy density 
given by 



H. 



x=l 



, s(a;) 



(32) 



, 45 Mp, ' ' ' 45 x3 
Y = n/s is the ratio of the total particle number density n to the entropy density s. The 
equilibrium number density Feq reads 



Y^,{x) = 4^^-^x'K2{x) . 
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(33) 



Pi = 1.0 




Pi = 0.1 




ai = 0.01 
ai = 0.1 
ai = 1.0 
ai = 2.0 


n,o = 6.4 X 10"^2 

°L 

n^o = 6.5 X 10-13 
n^o = 1.1 X 10-12 
n^o = 2.5 X 10-12 


ai = 0.01 
ai = 0.1 
ai = 1.0 
ai = 2.0 


n,o = 1.6 X 10-12 
n^o = 1.4 X 10-12 
n^o = 3.6 X 10-12 
n^o = 1.2 X 10-11 



TABLE IV: The relic number density n^o in terms of different ai and pi for the H2 case. 
In fact, 'Jeff is the reaction density defined by 

7e// = «eg (^e//^^) = -^e// (2) ^'l (^^ V^) , (34) 

with 

^eff = glff^E^E^aeffV^l-A/z , (35) 

where z = s/m^, Ki{x) and K2{x) are the modified Bessel functions. 

In our scenario, the exchanged particles may be Zi, Z2, h^, and if 2- is obvious 
that the case of exchanging gauge bosons Zi or Z2 is a p-wave annihilation process. If the 
exchanged particle is H^, the corresponding cross section will be suppressed by 1^^ /v\. For 
the case, the resonant condition 2m ^ myfi implies that the final states must be the Fermi 
pairs. In addition, the previous analysis indicates that the maximal cross section might be 
from the H2 exchanging process. Therefore, we study the and H2 cases in this subsection. 

Firstly we consider the H2 case. Due to the factor y^l — 4/z in Eq. (l35l) . we take 
m^^a/m? = 4.1 (namely e = 0.4878). At this point, 7e// becomes more larger than the 
rn^jjo/rn? = 4 case. Then we take pi = 1 (ps = 2.98). At this moment, 61 and 6^* may 
annihilate into h°h°, h^H^ and W2W2. Since and h^H^ have the similar form, the key 
quantity is a\ + ^2/8 for our calculation. Without loss of generality, one may take different 
values for ai and require ai = a2- The final results for different a\ have been shown in 
Table HVl In addition, we also calculate the p\ = 0.1 case, and list the corresponding results 
in Table HVl If p\ = 0.1, the final states have to be two SM Higgs bosons. In views of Table 
llVt we may find that the H2 case fails to suppress the relic number density of 61 and 61* . 

Now we assume the SM Higgs mass m/jO = 120 GeV in the case. Furthermore, one 
may obtain m = 59.3 GeV {epi = 3.5 x 10"^) from m^jjo/Tn? = 4.1. Because of m < 100 
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FIG. 4: Numerical illustration of the relic number density n^o as a function of 2m/m}^o near a 
resonance, where m^^o = 120 GeV has typically been taken. The dashed line denotes the present 
experimental upper bound on n^o . 



L 



GeV, the bound in Eq. (fT4|) is not valid. For m = 59.3 GeV, we take a'^^ < 4.6 x 10 ^^cm^ 
jiol l and derive the corresponding bound 



n^o < 1.1 X 10"^^ . (36) 



In this case, 6^ and 6^* mainly annihilate into the bottom quark pair. The annihilation 
cross section is given by 



/ . ^ ^ N 3 aim. s 



where is the bottom quark mass and F^jO ~ 3m/jom^/(87rfi;^) is the decay width of h^. One 
may obtain n^o = 1.2 x 10^^^ for ai = 2. This wonderful result indicates that our scenario 
may be consistent with the direct dark matter search bound. To illustrate, we plot the relic 
number density n^o versus the dark matter mass m in Fig. 4, where all annihilation channels 
have been considered. Using the results from CERN LEP-II, Datta and Raychaudhuri have 
derived m > 55.4 GeV 29|]. To show the resonance region, we choose 48 GeV < m < 72 
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GeV (0.8 < 2m/m/iO < 1.2) in Fig. 4. The peak around 2m/m/,o = 0.83 in Fig. 4 is 
due to the competition between and Zi resonances. For m/jo = 120 GeV, we find that 
56GeV ^ m ^ 60GeV can satisfy the requirement n^o < 1.1 x 10"^^. At this moment, 
one may obtain Qgo ^ 6.3 x 10"'', which is far less than the total dark matter density 



noAih^ = 0.111 ±0.006 



18|. 



VI. TWO HIGGS BI-DOUBLETS MODEL 

Motivated by the general two Higgs doublet model as a model for spontaneous CP viola- 
tion, one may simply extend the one Higgs bi-doublet LR model to a two Higgs bi-doublets 
LR model with spontaneous P and CP violation [l^ . Besides one left-handed Higgs triplet 
Al (3,1,2) and one right-handed Higgs triplet Ar (1,3,2), this model consists of two Higgs 
bi-doublets (2,2,0) and x (2,2,0), which can be written as 

The most general Yukawa interaction for quarks is given by 

-Cy = QI (n0 + Y,4> + F,x + F,x) Qr , (39) 

where Ql,b. = (""l,/?,, c^l,/?)^- Parity P symmetry requires Yg, Yg, Fg and Fg are hermitian 
matrices. When both P and CP are required to be broken down spontaneously, all the 
Yukawa couplings matrices are real symmetric. After the spontaneous symmetry breaking, 
two Higgs bi-doublets can have the following vacuum expectation values 

where Ki, K2, Wi and W2 are in general complex. Then we may obtain the following quark 
mass matrices 



M„ = —{YgKi+YgK2 + FgWi+FgW2] 



1 

Md = ^(^9«2 + YgKl + FgW2 + FgWl) . (AI) 

In the two Higgs bi-doublets model, the stringent constraints from the low energy phe- 
nomenology can be significantly relaxed. In Ref. |l3|, the authors calculate the constraints 
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from neural K meson mass difference Attik and demonstrate that a right-handed gauge 
boson W2 contribution in box-diagrams with mass around 600 GeV is allowed due to a can- 
celation caused by a light charged Higgs boson with a mass range 150 — 300 GeV. Therefore, 
we take vr ~ 2 TeV instead of the previous Vfi = 10 TeV for this section. It is worth- 
while to stress that our previous estimation is still right for this case except for the process 
^L^L* ~^ h^H^. {mao)fiOHO in Eq.(30) will be about 5 times larger than that in the vr = 10 
TeV case, which does not affect our conclusion. 

Since there are two Higgs bi-doublets, we can give more dark matter annihilation processes 
for 6l6T ^ HH and 6^61* ^ VH. In this model^ne may obtain three light neutral Higgs 
bosons and a pair of light charged Higgs bosons [30j]. The other Higgs bosons' masses are 
related to vpt- Although the annihilation cross section might be doubled or even increased 
by several times, it is still at least 10 times less than the direct dark matter search bound. 

An significant advantage of the two Higgs bi-doublets model is that the Yukawa couplings 
may become very large. In view of Eq. (41), one can explicitly understand this feature. For 
example, we require the couplings Yg and Yg are very large when wi ^ ki ^ k,2 ^ ^2- Then 
one may obtain more larger annihilation cross section for the 5^(5°* — // process than 
Eq.(23). For illustration, we take the maximal annihilation cross section for each quark pair 
final states 

„ 3 (aiYaWi)'^s 

^^-^^-^^ ~ s • 

where rrih denotes the mass of a light Higgs boson which comes from (p^. For at = 1, 
Wi ~ 246 GeV and m = 100 GeV, Yg ^ 4.3 can be obtained from Eq. fl20p when we take 
2m/mh = 0.8 and consider all quark final states but top quark. At this moment, we must 
consider the light Higgs h contribution to the direct dark matter detection experiments. The 
WIMP-nucleon cross section by exchanging h is given by 

^^^A^/^y (43) 

where / ~ Q.Q2YgWi/mu 15|- Using the above parameter setting, we may derive (j„ ^ 
6.1 X 10"^^ cm^, which is far more than the Zi exchanging case of Eq.(13). The more larger 
Yg is, the more larger cr„ is. Therefore, we can not give the desired relic number density 
through increasing the Yukawa couplings. 

Now we focus on the resonance case. For the H2 exchanging case, the results in Table HVl 
can be increased by about 5 times because of f^j ~ 2 TeV. On the other hand, more final 
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states would generally increase the partial width Tjjo_,fjij. Namely the case of more final 
states is equivalent to enhancing ai, which does no good for larger annihilation cross section 
as shown in Table IIV[ For the case, we may obtain the same conclusion as the one Higgs 
bi-doublet case. 



VII. SUMMARY AND COMMENTS 

In the Left-Right symmetric model with one Higgs bi-doublet, we have demonstrated 
that the cold dark matter constraints should be considered in a specific scenario in which 
the so-called VEV-seesaw problem can be naturally solved. In such a scenario, we find that 
(5° and are two degenerate and stable particles. To avoid the confiict with the direct 
dark matter detection experiments, we obtain the relic number density n^o < 4.8 x 10~^^, 
which implies that the two particles can't dominate all the dark matter. Subsequently, the 
lower bounds mcro > 0.13 GeV"^ and mcTo > 9.8 GeV ^ have been derived for the s-wave 
annihilation and the p-wave annihilation, respectively. In this paper, we examine whether 
our scenario can provide very large annihilation cross sections so as to give the desired 
relic abundance. We analyze in detail four classes of annihilation processes: 5^5^* — > //, 
^L^L* ~* ^l'^l* ~^ HH and 5^5°* VH . However, our analysis shows that this 

scenario fails to suppress the relic number density of 5° and 5£* except for the resonance 



case [35|J. For the resonance case, we obtain VL^o^h^ ^ 6.3 x 10~^, which is far less than 
the total dark matter density Qj^Mh'^ = 0.111 ± 0.006. Finally, we discuss the two Higgs 
bi-doublet model from the following three aspects: (1) ^ 2 TeV; (2) more final states; 
(3) large Yukawa couplings. It turns out that our previous conclusions can be generalized 
to the two Higgs bi-doublet model. 

In recent years, several authors have shown that it is far from natural for the minimal LR 
model to generate spontaneous CP violation with natural-sized Higgs potential parameters 
[sl, y, ^ . It is of importance for us to comment on some more general LR models with one 
Higgs bi-doublet [6|, 

mm 

, lll|. The differences mainly come from the complexity of 



Higgs potential parameter a2 and Yukawa couplings. We stress that our conclusion in Sec V 
could be generalized to these more general cases without any dramatic alternation because 
the gauge and Higgs sectors are basically the same. 
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APPENDIX A: SCALAR AND SCALAR-GAUGE TRILINEAR AND QUARTIC 
COUPLINGS 

We intend to calculate the cross section to the leading order for each process of dark 
matter annihilation. We first work in the framework of simple left-right symmetric model 
with one Higgs bi-doublet and one pair of LR triplets. To simplify our calculation, we take 
the decoupling limit in which n'^/vj^ ~ where = + 1^2^ denotes the EWSB scale. 
The VEVs of the Higgs bi-doublet are required to satisfy the low energy phenomenology 
constraint 1^2/ f^i < fnt/rrit, which may produce correct quark masses, small quark mixing 



angles and the suppression of flavor- changing neutral currents [31|, |32|, |33|, |3J] . For simplicity, 
we take K2 — which is a reasonable approximation at the leading order since is now 

around 10"^. Actually the limit K2 ^ brings additional advantage that the vacuum CP 
phase 62 could be taken zero safely without hampering the estimation. These approximations 
could largely simplify our calculation. 

The relevant scalar trilinear couplings and quartic couplings under the unitary gauge are 
shown in Table |Vl Here we write out the scalar-gauge interactions: 

^sisrvv = ^l^igWsL - g'Bf + gX^W^WE; (Al) 
^5ls,v = -w{Sld6l - 6ld6l)W^ + z6od6T{gW,L - g' B) + h.c.- (A2) 
Chvv = g'vn i^H'2[{gW3R - g'B){gWsn - g'B) + W^W^] + {-1=5],-W^W^ + /i.c.)| 

+ g^k\\^H^{W3RWt - W^lW^) - ]^{Hl + iAl)WEW^ + h.c. 
+ \h\W3L - Wm){W3L - Wm) + 2{WtWE + W^Wr )]^ (A3) 
Chhv = '-^H{dHt{W3L + W,R)-i5t+d5R~{gW-,R + g'B)-ig{5R-d5l-5ld5-R-)W+ 



T'g 
+ 1 

+ ^ 
2 



{A[dH{ - H:[dAl)W^ + h.c. 
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Interaction 


Coupling / vr 


Interaction 


Coupling / K 


Interaction 


Coupling 


Li Li Z 


P3 


Li Li 


ai 


Li Li 


ai 


Li L/ iX 




Li Li J- 




Li Li 1 


2a2 




ai 


/lO/lO/lO 


6A1 




«! + 03 




2a2 




6A4 


5l5l*A\A\ 

Li Li 1 -L 


ai + 03 




a\ + 03 




2A 


Li Lj 1 1 


ai 




a\ + 03 




6A4 




2P3 


Z 1 1 


Oi\ + Q3 


± ± 1 


2A4 


1^ L/ rt li. 


P3 




6pi 




2A' 








2(pi + 2p2) 


0-0 fj-O fj-O 


ai 
2a2 







TABLE V: The relevant trilinear and quartic scalar couplings, where the dimensional trilinear 
couplings with different scales vr and k are separated shown separately in two columns. A = 
Ai + 4A2 + 2A3 and A' = Ai - 4A2 + 2A3 



(A4) 



where the connection between weak eigenstates (W^sl, W^3_r, B) and physical states (Zi, Z2, A) 
are demonstrated by the following orthogonal transformation at the leading order: 

/ TA/.. \ / r„r n Qt,, \ ( 7..\ 

(A5) 



V B J 



cw sw 

-SiytiY C2W secvK Sw 
\ —\/C2wtw —tw \/C2W J 

The SU{2)l,r gauge coupling g and U{1)b^l coupUng g' are related to U{1)em gauge 
coupling e: 



Z2 



9 



9 



sin 9w vcos2^ 
Here our conventions are the same as those in Ref. 20|. 



(A6) 
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